Feasible optical weak measurements of complementary observables via a single 
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A general formalism for joint weak measurements of a pair of complementary observables is given. 
The standard process of optical three- wave mixing in a nonlinear crystal (such as in parametric down- 
conversion) is suitable for such tasks. To obtain the weak value of a variable A one performs weak 
measurements twice, with different initial states of the "meter" field. This seems to be a drawback, 
but as a compensation we get for free the weak value of a complementary variable B. The scheme 
is tunable and versatile: one has access to a continuous set of possible weak measurements of 
pair of observables. The scheme increases signal-to-noise ratio with respect to the case without 
postselection. 



The concept of weak measurements, introduced by 
Aharonov, Albert and Vaidman [![, gives us new insight 
and new tools to study quantum theory (for experimen- 
tal realizations see It allows us to have a fresh 
look at counterintuitive quantum phenomena [j| , and in- 
vestigate fundamental questions [fHEl . It is also a practi- 
cal tool for amplifying very tiny signals [lol - flij ]. difficult 
to observe in the conventional schemes. Very recently, 
with these ideas a direct measurement of the quantum 
wavefunction has been demonstrated 115], a protocol for 
amplification of weak charge signals in the context of a 
solid-state setup has been presented , and a scheme to 
amplify a cross-Kerr phase shift has been proposed [ItJ • 

Parametric down-conversion interaction (PDC), see 
e.g. [Til ], in the stimulated mode of operation, seems to 
be a good candidate for a process via which one could 
realize weak measurements. The coupling between the 
idler and signal beams is weak, and what is more impor- 
tant, it can be controlled by the strength of the pump 
beam, i.e., can be made as weak as we wish. However, 
the interaction Hamiltonian is at first glance totally un- 
suitable for this purpose, as it does not have the usual 
form of a coupling between a system observable (A) with 
an observable (p) of the "measuring device" (or rather 
measuring system) , such as g A <£> p, where g denotes the 
coupling strength. In this letter we propose the follow- 
ing. One can treat the idler mode of the parametric pro- 
cess, denoted by ad, as the "device" or "meter" subsys- 
tem, and the signal mode as the system to undergo a 
weak measurement (a s ). If one as usual describes the 
(pulsed) pumping fields as classical, which is a standard, 
well working, approximation for narrow crystals and a 
not too strong pump, the parametric down conversion 
Hamiltonian can be put as H = gf(t — to)(a s ad + aja^), 
where / represents a temporal overlap function of the 
pulse with the volume of the crystal (we can safely as- 
sume that it is very close to a Dirac's delta 5(t — to)), and 



g is an effective coupling defined by the pump strength 
and geometry, second order non-linearit y an d crystal size 
(for a derivation and details see e.g. [la]). This form 
of interaction does not have the characteristic form re- 
quired for weak measurements. However, if one puts 



B = (a s + aJ)/-\/2 and A = i(a s — a\)/y2, and for the 

a d )/V2, the 



"meter" q = (a d + a d )/y/2 and p = i(a\ 
Hamiltonian acquires a new form 



U = gS{t - i )(A ® p + B ® q), 



(1) 



i.e., as the sum of two Hamiltonians used for weak mea- 
surements, which apply to two complementary situa- 
tions, both for the system and the "meter". Notice 
that [q,p] = i, and [A, B] = i. We shall show be- 
low that, surprisingly, a Hamiltonian of such a generic 
structure still allows weak measurements, and what is 
promising, it is very versatile: one can also put it as 
% = gS(t — to) ( A' ® p' + B' ® q') . The new operators are 
suitable pairs of quadratures B' = (e~ %v a s + e %v a\)/^/2, 
A' = i(e-^a s - e^a\)/^, q' = (e^a d + e~ il f a\) / y/2 
and p' = i(e~ lip a d — e l<p ad)/ V%, for a given phase ip. We 
have the same structure for every tp\ We shall always put 
ip = 0, however all that will be said below also equally 
pertains to other possible values. 

Our claim is that a PDC-type Hamiltonian allows weak 
measurements of pairs of (complementary) observables. 
Depending on the pairs of observables (defined by ip), 
one has to prepare different initial states of the " measur- 
ing" subsystem, and different selection procedures. Still, 
the coupling remains the same, thus we have in our labs 
a ready, versatile, tool for weak measurements. To ob- 
tain the weak value of a variable A one has to perform 
weak measurements twice (which is neither conceptually 
nor practically difficult). This seems to be a drawback, 
but as a compensation we get for free the weak value 
of the complementary variable B. Furthermore, a weak 
measurement of a pair of observables may lead to a sig- 
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nificant amplification of signal-to-noise ratio (SNR) with 
respect to the case without postselection (provided quan- 
tum noise, or shot noise, is dominant). This is in sharp 
contrast with the case of weak measurements of single 
observables. For such cases, once quantum noise is dom- 
inant and the probability decrease due to postselection 
is taken into account, SNR and measurement sensitivity 
cannot be effectively improved with respect to the case 
without postselection [19(. 

Let us remark that a Hamiltonian like ([T]) cannot gen- 
erally lead to a decoherence of the state of the system 
into an eigenbasis of neither A or B. This may be the 
case, e.g., for initial Gaussian states of the meter system. 
From the point of view of measurement theory it does 
not provide an eigenselection, and thus is useless - we 
stumble here on a new feature distinguishing weak and 
standard measurements. 

Let us present a general description of a weak mea- 
surement of two noncommuting observables of a system, 
A and B. We assume that the interaction Hamiltonian 
of the system s and a measuring, or "meter" , system d 
is given, as above, by (|T|). In initial steps of the follow- 
ing discussion neither the pair (A,B) nor the pair (q,p) 
needs to be a complementary pair of observables. We put 
h = 1 and assume that all factors in gA <E> p and gB ® q, 
especially g, are dimensionless. We assume throughout 
that the adjustable constant g, the effective strength of 
the system-meter interaction, is small enough, so that 
the weak interaction conditions [20j <?||A||Af> <C 1 and 
g||Bj|A(7 <C 1 are satisfied. Thus, in perturbation expan- 
sion we keep only terms of the lowest order in g. There- 
fore, the time-evolution is given by 



JJ _ e -ig(A®p+B<»g) 



1 



ig(A®p + B®q). (2) 

Let |i/>j) and |0) denote the initial states of the system 
and of the meter, respectively. After the interaction, the 
state of the compound system is given by 

U \4>) « \1>i) \4>) ~ *5(A ® p + B ® q) \rPi) |0) . (3) 

The subsequent postselection is defined by projection 
onto the state of the system \ipf). Upon a successful 
postselection, the final state (not normalized) of the mea- 
suring system is given by 

|^) = ^/|C/|^}|^) 

« (^Mi) |0) - ig((ipf \ A \i/ )i )p+ (^| B |Vi) q) |0) 
= (fl>f\il>i){l-ig(A wP + B m q)}\4>), (4) 

where A. w and B w are the weak values defined according 
to 



(#|A|^) 



(5) 



(ipf\ipi) 

and similarly for B w . The success probability of the post- 
selection is given by 

P={4>'W) 

« \{ipf\tpi}\ 2 (1 + 2g (p) JmA w + 2g (q) ImB w ) . (6) 



The average shift of the pointer variable q, conditional 
on preselection to \ipi) and postselection to \i/jf), is given 
by 



Sq = 



WW) 

2glmR w YBxq + gRcA w (—i[q,p\) 
+glmA w ({q-(q),p-(p)}), (7) 



where [t, s] and {t, s} denote the commutator and anti- 
commutator of t and s, respectively, (s) denotes the ex- 
pectation value of s in the initial state, van? = (q 2 ) — (q) 
denotes the variance of q, or (Aq) 2 , in the initial state. 
The average shift of the pointer variable p is similarly 
given by 



Sp - 



2gImA„,varp — gRcH w (—i[q,p]' 
+gImB w ({q- (q) ,p- (p)}) . 



(8) 



The approximations in calculating average pointer shifts 
are legitimate if the higher-order terms are negligible, 
therefore we further require that <?||A||Ap <C | (ipf\tpi) \ 
and p||B||Ag <C \{ipf\ipi)\. Obviously, preselected and 
postselected states should not be orthogonal. 

Let us now assume q and p are a pair of conjugate vari- 
ables of the measuring device, i.e., [q,p] = i. Therefore 

Sq 2gImB,„var(7 + gReA w 

+glmA w ({q-(q),p-(p)}) (9) 
Sp ~ 2gImA,„varp — gRcB w 

+glmB w ({q-(q),p-(p)}). (10) 

For simplicity, we additionally assume that the meter is 
initially in a Gaussian state 



<pAq(q) 



1 



(V^Aq) 1 / 2 



Thus, the anti-commutator terms in Eqs. (|9I10I) vanish. 
Therefore 



Sq ss 2glmB w (Aq) 2 + c/ReA„ 
Sp « 2glmA w (Ap) 2 - gRcB^ 



(12) 
(13) 



Note that, Ap = l/(2Ag) for the Gaussian state in (jlll) . 

With the same interaction Hamiltonian but two differ- 
ent initial states (f)Aq 1 {q) and 4>Aq 2 ( c l)i we have 



Sq! = 2gImB w (A qi ) 2 + gReA v 
Sq 2 = 2gImB w (Aq 2 ) 2 + .gReA„ 
Spi = 2,gImA Ii ,(Api) 2 - gReB t 



5p 2 = 2glmA w (Ap 2 ) 2 - gReB % 



(14) 
(15) 
(16) 
(17) 
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Thus, 



2gImA„ 

gReA v 
gRcB v 



6pi - Sp 2 



(Api) 2 
Sqi 



- (AP2) 2 
Sq 2 



(Agi)a - (Ag 2 ) 2 
(A qi ) 2 6q 2 - (Aq 2 ) 2 S qi 

(A gi ) 2 - (Afe) 2 
(A P2 ) 2 S Pl - (A Pl ) 2 Sp 2 

{ApiY - (A P2 ) 2 ■ 



(18) 
(19) 
(20) 
(21) 



Hence, with a pair of weak measurements, each with a 
different initial state of the meter, however both using the 
same Hamiltonian, or if one likes the same experimental 
setup, one can determine the real and imaginary parts of 
the weak values of a pair of non- commuting observables. 
The two initial states of the meter system must differ in 
spreads of pointer variables p and q. 

In real experiments, postselection may not be a pro- 
jection onto a pure state, and the preselection might not 
be ideal. We extend the discussion to such cases in the 
Supplemental Material. 

Thus far we have discussed the unexpected features 
of a PDC Hamiltonian. Let us now propose one of 
the possible blueprints for an operational realization of 
such a scheme, employing standard tools of modern 
quantum optics, see Figure 1. With the laser pump 
beam on, the interaction within the crystal is described 
by (fT]). Our "system" and its initial state can be a 
weak coherent state of light, |a). It is incident on a 
beam splitter which is slightly unbalanced. The out- 
put light from the beam splitter is described by \ipi) = 

-9L(1 - e )) + e)) , where < e < 1. This will 



be our initial state of the system. The exit beam splitter 
in the system interferometer is a symmetric 50-50 one. 
The meter system, or "pointer" beam d, is fed with a co- 
herent or squeezed state. The proposed postselection of 
the system states corresponds to a detection of a single 
photon at the port 3s (and no photon at port 4s). The 
interferometer is tuned such that if the pump beam is off, 
and thus there is no interaction within the crystal, and 
e = 0, only detector 4s registers light, i.e., 3s is a dark 
output in such a case. Formally this means that we se- 
lect states \tf, f ) = |l) 3s |0) 4s = ^(|l) s , |0) s - 1 |0} s , |1)J. 

In such a case the weak values of A = -j^(a s — a J) and 



B 



1 

V2 



(/Is 



4) 

A ( , 
B„ 



are given by 



V2 K 



1 a 

~2ea ~ 2 ' 

i ia 

2ea + T 



1 

2ea' 
i 

~2sa' 



(22) 
(23) 



Equations (|T2"j) and (fi"3"|) give the average shifts of the 
"pointer" beam in terms of q and p. For the presented 



scheme they read: 
5q = 

Sp - 



g Kea 
e \a\ 2 



(Aqf 



g Ima 



(Ap) 2 



(24) 
(25) 



The shifts Sq and Sp can be measured with a homodyne 
detection. Parameters Aq and Ap depend on the ini- 
tial state of the pointer beam. As Sq and Sp are di- 
rectly accessible, and e is also known and controllable, 
this scheme gives, for example, a direct method to mea- 
sure the (small) coupling coefficient g. 

Generally under suitable conditions one can achieve 
with weak measurements an amplification of signal and 
signal-to-noise rate. If the final state of the system is 
not postselected, the average shifts of the pointer vari- 
ables q and p are Sqo « g and Spo ss g (B) i , where 
(fi) i = (ipi\Q\ijji) denotes the expectation value of f2 in 
the initial state \ipi), which always lies inside the range of 
the eigenspectrum of f2. By comparing these values with 
([121131) . we can calculate the signal amplification factors 
for our weak measurement scheme, defined as K„ = -P- 

' q Sqo 



and Kn 



As we shall show, each of them can be 



Sp 

lp = 

much larger than 1. The reason for this is that the weak 
values could lie far outside of the eigenspectrum, and Aq 
of the initial state is an adjustable parameter. 

If no postselection is employed in reading of the shift 
of the pointer variable q, the signal-to-noise ratio (SNR) 



is given by lZ q = 



where N denotes the number of 



q Aq/VN' 

repeated readouts, whereas the SNR for weak measure- 
ments is given by H' a = A J q /ATD , where P « | W>j#») 1 2 

To 



1 Aq/VNP 

stands for the success probability of postselection 



calculate 1Z' q we use the initial value of Aq. This is jus- 
tifiable as in the lowest order of perturbation series, the 
standard deviation of q of the pointer is not changed by 
the postselection. Therefore, the amplification factor of 

the SNR A q = ^ i s 



2ImB tu (A<7) 2 + ReA-u 



MM) 



(26) 



In the case of the shift of p, due to the postselection in 
the weak measurement, the SNR is amplified by a factor 



2ImA u ,(Ap) 2 - ReB„ 



(27) 



To illustrate the above, and to show capabilities of the 
proposed experimental scheme, let us calculate in detail 
the amplification factors for the setup of fig. 1. Without 
the postselection the average shifts of the pointer position 
and the pointer momentum, as measured via a homodyne 
detection, are given by Sqo ~ — gRea and Spo w glma. 
If the initial state of the pointer is a coherent one, we 
have (Ag) 2 = (Ap) 2 = ±, whereas Sq = -f raf and Sp = 



4 




FIG. 1. (Color online) A sketch of a possible experimental 
setup. On a non-linear, say, BBO crystal impinges a pump 
beam of a power tuned such that the three- wave mixing inter- 
action leads effectively to the coupling Hamiltonian |T]) with 
a small enough g. Path d represents the "pointer" beam (de- 
vice), whereas paths s and s' represent the measured system 
(beams are phase matched with the pump beams, which in 
turn is treated as classical). The exit beamsplitters are sym- 
metric 50-50 ones. The upper input beamsplitter is slightly 
biased (see the text), whereas the lower one should be optimal 
for the homodyne measurements. In the path d leading to the 
BBO crystal one can put additional, active or passive, optical 
elements modifying the properties of the "pointer" beam. 

- Is^f . All this shows that in the weak measurements the 

6 \a\ 2 

measured shifts are amplified by a factor of K = jj^p- 
This amplification factor could be very large for a small 
£ |a| 2 . 

Let us move to the signal-to-noise ratios. Using ([26l 
l27f one can establish that the SNRs are amplified by a 

factor of A = Ky/P = e A 1 , which approximately 
reads A s» tK for |a| -C 1. Therefore, by choosing a 
very small e, we can have a very large signal amplifi- 
cation, while surprisingly by making the system beam 
very weak (i.e., a very small |a|), both the shift and 
the SNR in the meter beam can be amplified. How- 
ever, the amplification cannot be infinitely large. For 
the high-order terms of g to be negligible, and thus all 
the approximations legitimate, one must require that 
.9 ^ I (V'/IV'*) I ~ £ \ a \- The overall success probability 
of postselection is P m \ (ipf\ipi) | 2 ~ e 2 |a| 2 . We use 
small a to achieve large SNR. This is especially useful in 
cases in which the input beam is technically limited to 
very low number of photons, or the quantum noise needs 
to be suppressed in order to register observable pointer 
shifts. 

For a squeezed state of the meter, the shifts Sq or Sp 
will be further amplified. From eqs. (I24l25[) . for Aq 3> 
1/2 we have Sq « -f ^(Aq) 2 , and for Ap > 1/2 we 

have Sp m ^jSp-(Ap) 2 . For the aforementioned Aq the 
signal Sq is amplified by a factor of K q = , and 

SNRby A q = K q VPn£f. 

A typical thin BBO crystal, and a moderate pump 
power, lead to g around 10 -6 . If the meter beam is ini- 
tially in a coherent state, for e ~ 10~ 2 and a ~ 10~ 2 , we 
have an amplification of the signal by 6 orders of magni- 



tude and an increase of SNR by 2 orders of magnitude. 
If the meter beam is prepared in a squeezed state, with 
Ap or Aq squeezed by a factor of 10, a further ampli- 
fication of the signal and SNR by additional two orders 
of magnitude is possible. Other sources of noise, e.g., 
systematic errors linked with alignment in the interfer- 
ometers, etc, may reduce the actual amplification factor. 
Note that the condition g <C | (i/jf\tpi) | ~ e\a\ is satisfied 
for thin crystals with a not too strong pump. In such a 
case, corrections to higher orders of g is not necessary. 

Finally, let us remark that a scheme for a weak mea- 
surement involving a PDC Hamiltonian was proposed in 
[2lj |. However, in this case, the coupling was strong, the 
set-up worked in an optical parametric amplifier regime, 
and weak measurements were due to a weak coupling of 
the detector to the field in the idler mode. Thus the inter- 
action (Q} was not directly responsible for the weak mea- 
surement, but rather used for preparation of initial states. 
What we propose is absolutely different. The PDC in- 
teraction Hamiltonian, in the weak coupling regime (a 
very small g), which at first glance looks unsuitable for 
an interaction, leads to weak measurement that has un- 
expected properties and is highly versatile. It allows a 
quantum optical scheme for weak measurement of a pair 
of complementary observables. We can also choose which 
pair we are interested in. The schemes described above 
are feasible, as the current methods of handling the PDC 
process are very well developed. The general scheme for 
weak measurement of a pair of observables can also be 
implemented by other processes endowed with a simi- 
lar interaction Hamiltonian, such as the opto-mechanical 
interaction in the strong coupling regime, and the four- 
wave mixing process in cold atoms. 
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SUPPLEMENTAL MATERIAL: A DISCUSSION 
ON IMPERFECT PRE- AND POSTSELECTION 



General case 

In real experiments, the postselection may not be a 
projection onto a pure state of the system, e.g., detec- 
tors not sensitive to photon numbers can be modeled as 
a projection onto the subspace with at least one photon. 
Suppose the initial state of the system is prepared in a 
general mixed state p s and the postselection is a gen- 
eral projection onto a subspace, denoted by II/. With 
the evolution governed by (1) in the main text, the 
state of the meter after the interaction and the posts- 
election is given by p' d — tr s [UfU(p s ® p d )U^], where 
U is the time evolution operator due to the interac- 
tion, p d is the initial state of the measuring device, and 
the trace is taken over the Hilbert space of the system. 
The average shift of a pointer variable M conditional 
on the pre- and postselection of the system is given by 
5M = tr d {Mp' d )/tr d {p' d ) - tr d (Mp d ). After a detailed 
calculation, we get exactly the same formulas (7-8) as 
for the case of pure pre- and postselected states in the 
main text. Of course, the weak values (A„, and ~B W ) for 
general preselection p s and general postselection 11/ are 
defined according to [2(| 



_ tr(U f n Pa ) 
{ )w ~ tr(U fPs ) • 



where fl = A or B. This is a natural extension of the 
original definition of weak values. A further generaliza- 
tion to positive operator-valued measurements (POVM) 
is possible, see [221 ] . however we shall not address it here. 



Our example 

Let us consider how the results in our proposed ex- 
periment can be affected, if the post-selected state is 
not a pure one. Photons from the dark port 3s might 
be collected by a photon detector which cannot distin- 
guish photon numbers, which is a quite typical situation. 
The postselection is therefore better modeled as a pro- 
jection onto a subspace with at least one photon, i.e., 
IT/ = 7 3s — |0) 3s (0|. The pre-selected state can be writ- 
ten as = ^=(1 - e) 

beamsplitter this leads to |— ea) 3s \ia) As 
variables of ports 3s and 4s, we have 



^|(1 + e)j .At the output 
In terms of the 



A = Tl^ 5 ~ °^ = 2 a3s + a ^ + ^ a4s ~ a ^ 



and 



B = + al) = ^ -i{a 3s - 4j + a 4s + a\ 



Thus, from (1281) we obtain the weak values 



4 s 



B„ 



1 

2ea 
i 

2ea 



Rea 



— Ima 



1 

'~2ea' 
i 

~2sa 



(29) 
(30) 



(28) 



which barely differ from Eqs. (22) and (23) in the main 
text, and thus can be approximated to the same values. 
The correction is negligible for e -C 1, therefore the av- 
erage shifts of q and p are still given by (24) and (25) in 
the main text. 

For an even better modeling of the photon detection, 
one may take into account the cases in which the detec- 
tor is not clicking, while there is one or more photons. 
Therefore, one can write 11/ = 77 (I 3s — \0) 3s (0|), where 
77 is the detector efficiency. Even with the detection effi- 
ciency taken into account, all the above results still hold 
since we only consider the data from a successful posts- 
election. The presence of an inefficiency only decreases 
the probability of a successful postselection. 



